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We propose a class of observables constructed from lepton energy distribution, which are inde- 
pendent of the velocity of the parent particle if it is scalar or unpolarized. These observables may 
be used to measure properties of various particles in the LHC experiments. We demonstrate their 
usage in a determination of the Higgs boson mass. 
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The data and analyses from the experiments at the 
CERN Large Hadron CoUider (LHC) are attracting in- 
creasing attention. The main goals of these experiments 
are discoveries of the Higgs boson and of signals of physics 
beyond the Standard Model. Once the Higgs boson or 
other new particles are found, the next step is to uncover 
properties of these particles. There are well-known chal- 
lenges in performing accurate measurements for such pur- 
poses: (i) In reconstructing the kinematics of events, it is 
difficult to reconstruct jet energy scales accurately. This 
is in contrast to electron and muon energy-momenta, 
which can be measured fairly accurately, (ii) Often inter- 
esting events include undetected particles which carry off 
missing momenta. In this case, reconstruction of miss- 
ing momenta is non-trivial, (iii) The parton distribu- 
tion function (PDF) of proton is needed in predicting 
cross sections and kinematical distributions. Our cur- 
rent knowledge of PDF is limited, which induces rela- 
tively large uncertainties in these predictions. As a con- 
sequence of these difficulties, for instance, it is difficult to 
measure the energy-momentum of a produced new par- 
ticle or predict accurately their statistical distributions. 

To circumvent these difficulties, many sophisticated 
methods for kinematical reconstruction of events have 
been devised. See, for instance, [l|, Q and references 
therein. In these methods, one takes advantage of var- 
ious kinematical constraints, which follow from specific 
event topologies, to complement uncertainties induced by 
the above difficulties. Still, in most cases challenges re- 
main to reduce systematic uncertainties originating from 
ambiguities in jet energy scales, PDF and other non- 
perturbative or higher-order QCD effects. Generally it is 
quite non-trivial to keep these uncertainties under con- 
trol, both theoretically and experimentally. 

In this paper we propose a class of observables which 
can be used to measure properties of particles produced 
in the LHC experiments, by largely avoiding the above 
uncertainties. We consider a particle X which is scalar 
or unpolarized and whose decay daughters include one 
or more charged leptons i= e or n . The observables 
are constructed from the energy distribution of £ and are 
independent of the velocity of X. Although experimental 
cuts and backgrounds induce corrections to this property, 
we will show that systematic uncertainties are suppressed 
and can be kept under control. Furthermore, we can 
utilize the degree of freedom of the observables to reduce 
or control effects induced by cuts and backgrounds. 



In the first part of the paper, we explain the construc- 
tion of these observables, in the two-body decay and 
multi-body decay cases separately. In the latter part, 
we demonstrate its usefulness in a determination of the 
Higgs boson mass using the vector-boson fusion process. 

2-body decay: X ^t + Y 

Suppose a parent particle X decays into two particles, 
one of which is a lepton i whose energy is Eq (monochro- 
matic) in the rest frame of X. In the case that X is 
scalar or unpolarized, the normalized lepton energy dis- 
tribution in a boosted frame, in which X has a velocity 
/?, is given by 



Vp{Ee;Eo) = 



1 



{ev - e~y)Eo 



'^Ea <Ei< e'^So), (1) 



in the limit where the mass of the lepton is neglected. 
Here, y is the rapidity of X in the boost direction [13|. 
related to /3 as e'^y = {1 + f3)/{l - [3). The step function 
is defined such that 6{cond.) = 1 ii cond. is satisfied, and 
9{cond.) = otherwise. 

We construct an observable g{Ee/Eo) from the lepton 
energy Ei in the boosted frame such that the expectation 
value (g) is independent of /?. For later convenience, we 
write g{Et/Eo) = [dG{x)/dx]x=Ee/Eo- Hence, 



{g)= dEeVp{Ee;Eo)g{Ei/Eo) 



G{ey)-G{e'y) 



(2) 



This shows that {g) is the even part F{y) + F{—y) of 
F[y) = G{ey)/{ey - e-y). It is ;3(y)-independent when 
F is an odd function of y plus a constant independent of 
y. Namely, G{ey) = (even fn. of y) + const, x (e^ — e^^). 
Let us demonstrate usage of (g). In the case G = 
ey - e-y, we obtain (l/E^) = l/E^ Thus, (mathe- 
matically) we can reconstruct Eq from the lepton energy 
distribution irrespective of /3 of the parent particle. In 
the case G(ey) = (even fn. of y), {g{Ei/Eo)) = 0. Con- 
versely, we can adjust Eq, which enters g{E£/Eo) as an 
external parameter, such that (g) vanishes and determine 
the true value of Eq. We give two examples of G for the 
latter case: (a) Gn — l/[2ncosh{ny)], corresponding 
to g = E^+^E^-^Ei"" - S2")/(£;2» + S|")2. As n in- 
creases, contributions of the Ei distribution from large \y\ 
region {Ee ^ Eq or Ee ::!> Eq) become more suppressed. 



(b)G, 



(fc) 



(r < e 



-|y| 



< l)(e- 



-\y\ 



r), corresponding to 



g = 0{r < Ee/Eo < 1) - {Eo/Ee)^9{l < E^/Eo < r'^). 



This observable is independent of the Eg distribution in 
the regions Ei < tEq and Ei > r^^Eo. 

Many-body decay: X -^ i + anything 

We assume that we know the theoretical prediction for 
the lepton energy distribution in the rest frame of the 
parent particle X, dT x^i+a.nything/ dEo , where Eq is the 
lepton energy in this frame. We define an observable in 
the boosted frame as 



Og=M 



dEp drx^^+anything 



El 



dEo 



dx 



G(x) 



.(3) 



x=Et/Eo 



{Og)d = to determine mx- Generally various experi- 
mental cuts, which affect D{Ee), are imposed due to de- 
tector acceptance effects and for event selection purposes. 
Furthermore, there are contributions from background 
processes which also modify D{Eg). After incorporating 
these effects, there is no guarantee for (Og) to vanish. 
Suppose that D{Ei) is modified to D{Ee) + 5D{Ef) by 
these effects, where the distributions are normalized as 
JdEeD ^ JdEe{D + SD) = 1. In the case \5D/D\ < 1, 
if {Og)d+sd = is used to extract mx, the obtained 
value is systematically shifted from the true value by an 
amount 



It depends on the lepton energy Ef^ in the boosted frame 
and on the parameters of dT/dEg such as the parent par- 
ticle mass mx ■ A/" is an arbitrary normalization constant 
independent of Ei. We can prove (see below) that, if X 
is scalar or unpolarized, and if we take the same G as in 
the 2-body decay case, (Og) is independent of f3 of X. 
In particular, in the case that G(e^) is an even function 
of y, (Og) — 0. The parent particle mass enters as an 
external parameter in the definition of Oq- Hence, we 
can use this property to determine m,x, provided that 
other parameters are known, g in the 2-body decay case 
can be regarded as a special case of Og- 

Proof: The lepton energy distribution in the boosted 
frame is given by 



MEi) - 
ce. 


jdE','- 


i X->f-|- anything 

dEl, 



-Dp^E.-E'^) 



(4) 



(Og)= / dEifp{Ei)OG 

^, dr dr G{eyE'jEo)-Gie'yE'jEo) 



= U / dEpdEQ 



dEo dE'^ 



Eq E'o {ev - e-y) 



(5) 



where we integrated over Ef^. In the case G{x ^) = 
G(x), GieyE'jEo) - Gie-yE'jEo) = G{eyE'jEa) - 
G{eyEo/EQ) is anti-symmetric under the exchange of Eq 
and Eq, while all other parts are symmetric. It fol- 
lows that (Og) = 0. In the case G{x) — x ~ a;^^, 
GieyE'jEo) - Gie-yE'jEo) = {E'JEq -f Eo/E^){ey - 
e ^), so that y-dependence cancels out. (Q.E.D.) 

Since we must know the theoretical prediction for 
dr/dEo to construct Og, we consider use of Og mainly 
for the purpose of precision measurements, after the na- 
ture of the parent particle (such as its decay modes) is 
roughly determined by other methods. In principle, dif- 
ferent functions G in constructing Og may be used to 
determine simultaneously more than one parameters in 
a decay process. In this first study, however, we consider 
the case where only one parameter is unknown. 

Under a realistic experimental condition, in which /3 
of X has a distribution Dx{/3), the lepton energy dis- 
tribution in the laboratory frame is given by D{Ef) = 
J dl3 Dx{P)fj3{Eg). In this case, the expectation value 
{Og)d = JdEeD{Ee)OG is independent of Dx- Here- 
after, we choose G to be an even function of y and use 



5mx = -{Og)sd 



'/OOgX 
\dmx / , 



(6) 



where we neglected 0{SD^) corrections. We can use this 
formula to study systematic corrections in the determi- 
nation of mx ■ SD should be estimated by Monte Carlo 
(MC) simulations which take into account realistic ex- 
perimental conditions. Errors in the estimate of SD con- 
tribute _as_ systematic uncertainties in the determination 
of mx- 



IM 



The production of Higgs bosons via vector-boson fu- 
sions is expected to be observed with a good signal-to- 
noise ratio, if the Higgs boson mass m,H is within the 
range 135 GeV ^ mn ^ 190 GeV, hence it is considered 
as a promising channel for the Higgs boson discovery. 
We perform a MC simulation to study feasibility of m.H 
determination using the observable Og and the decay 
modes H -)■ VFVF^*^ -^ Uvv {U = e^,/i/i,ee). We gen- 
erate the events for the signal and background processes 
using MadEvent Q, which are passed to PYTHIA [1] 
and then to the fast detector simulator PGS jSj. We set 
Vi = 14 TeV. 

The strategy of our analysis follows, to a large extent, 
that of p, which studied the prospect of Higgs boson 
search using the vector-boson fusion process in the AT- 
LAS experiment. We first repeated the analysis of [y] in 
the case oi H ^ e^xvv mode using our analysis tools and 
imposing the same cuts. We reproduced the numbers in 
Tab. 4 of Q reasonably well, considering differences in 
both analyses, such as different detector simulators, dif- 
ferent PDF's, and different jet clustering algorithms (we 
use the cone algorithm with R — 0.5; we do not correct 
the jet energy scales given by the output of PGS) : regard- 
ing the signal events, we reproduced the efficiencies of the 
cuts involving only leptons within a few % accuracy and 
those involving jets within a few tens % accuracy; in to- 
tal, the efficiency of all the cuts was reproduced with 6% 
accuracy. On the other hand, the cross section of the sig- 
nal is smaller by 30% in our analysis as compared to that 
of [g]. The difference originates from the different scales 
and different PDF's used in the event generators. We do 
not correct the difference in normalization of the cross 
sections; it may result in overestimates of the statistical 
errors given below. 

For simplicity, in our analysis we omit the Higgs boson 
production via gg fusion and the Higgs decay modes in- 
cluding t's. In principle, part of these modes including e 
or /i in the final states can be used as signal events, since 
dF/f^^+anything/rf^o is Calculable. 



efj, mode [fb] 
Hfi mode [fb] 
ee mode [fb] 



Signal (mu) 
150 GeV 180 GeV 200 GeV 



2.00 
1.33 
0.76 



2.21 
1.47 
0.89 



1.07 
0.68 
0.42 



Background 

tt + Wt WW +]ets 



0.53 
0.31 
0.24 



0.13 
0.05 
0.04 



TABLE I: Cross section x efficiency, after all the cuts. e,/i 
from r are not included. 



In our analysis of mn reconstruction, ideally two crite- 
ria need to be satisfied to ensure use of the observable Oq'- 

(1) The lepton energy distribution in the Higgs rest frame 
agrees with the theoretical prediction (ITh-^Uvv/ dE^. 

(2) The lepton angular distribution in the Higgs rest 
frame is isotropic. The effects of cuts should not violate 
these criteria significantly. 

Jets in the signal process are associated with the Higgs 
production process and are independent of the Higgs de- 
cay process. Therefore, cuts involving only jets would 
not affect the above criteria but only affect the (5 distri- 
bution of the Higgs boson. By contrast, cuts involving 
leptons can affect the above criteria significantly. 

Taking this into account, we use the same cuts as in 
except the following two cuts: 

(A) Lepton acceptance 

No muon isolation requirement; 

Pxie) > 15 GeV, Prifi) > 10 GeV, |r;f| < 2.5; 

Only one pair of identified leptons in each event. 

(B) Lepton cuts 

Mu < 45 GcV, Prie), Priti) < 120 GeV. 

Concerning (A), since lepton isolation requirement as 
well as cuts on lepton Pt and rj bias the lepton angu- 
lar distribution in the Higgs rest frame, we loosen the 
cuts and requirement as much as possible. Concern- 
ing (B), it is important to select events with leptons 
in the same directions, in order to reduce background 
events. We note that, if we impose the same cut on 
Mil in the theoretical prediction for dTH^Uui^/dEo, this 
cut does not contribute to systematic errors. Hence, we 
tighten the Lorentz invariant cut and omit other frame- 
dependent angular cuts. By modification of the cuts (A) 
and (B), the efficiency of the signal events reduces by 
a few tens %, while variations of the efficiencies of the 
background events are small. (l5| 

Tab.|I] lists estimates of (cross section) x (efficiency) of 
the signal and background events, after all the cuts 
are imposed. |16j For the backgrounds, we simulate only 
ti + Wt and M^W-l-jet (electroweak) events, which are 
shown to be the major backgrounds in [6|. We estimate 
contributions of other backgrounds to be negligible, com- 
pared to the uncertainties discussed below. Consistently 
with the signal, we omit background events with £ ~ e, ii 
from T. 

We test the two criteria with the signal events which 
passed all the cuts. The histogram in Fig. 1 shows the 
lepton {£ = ij,) energy distribution in the Higgs rest frame 
for the H -> ^^i/v mode and IVIC input m|J^ — 150 GeV, 
where we looked up the parton-level neutrino momenta 
in each event to reconstruct the Higgs momentum. We 
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FIG. 1: Lepton energy distribution in the Higgs rest frame. 
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FIG. 2: Lepton cos 9 distribution in the Higgs rest frame. 



generated a large-statistics event sample (about 12,000 
events after applying all the cuts) in order to focus on 
the systematic effects caused by the cuts. The theoretical 
prediction for dV h-^Uw / dE^ with the cut Mu < 45 GeV 
(no other cuts are applied) is also plotted with a red line. 
A good agreement is observed. 

Figs. 2 show the lepton cos6' distributions in the Higgs 
rest frame of these events, where is the angle measured 
from the Higgs boost direction. The total event sample 
is divided into four groups of equal size, in the increas- 
ing order of /? (or y) of the Higgs boson. The average y 
value for each group is displayed. The events with low 
y have a distribution closer to isotropic one, whereas the 
events with large boost factors have strong distortion of 
the cos 9 distribution. Since Higgs bosons with large y 
are boosted mostly in the beam directions, the lepton 
Pt and ij cuts and electron isolation requirement bias 
the cos 9 distribution. In particular, depletion of events 
in the cos ~ 1 region of large y samples is caused by 
the lepton rj cuts. On the other hand, Higgs bosons with 
small y are boosted in random directions, so that the 
lepton acceptance effects do not bias the lepton cos 6 dis- 
tribution strongly. Thus, the second criterion is satisfied 
only by events with small boost factors. 

We may choose G appropriately to suppress contribu- 
tions of events with large y. By examining various G, we 
find that G = G„ for n « 4 is an optimal choice. As 
n increases, (Og) becomes less sensitive to events with 
large y but more sensitive to statistical fiuctuations.flTj 

In niH determination, typically there are two solutions 
which satisfy {Oq} = 0, one above and one below the 
WW threshold 2Mw- We befieve that, unless niH is very 
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FIG. 3: Muon energy distribution after applying all the cuts 
for H — >■ ^^Li^v mode. Histograms represent muons from 
the signal (corresponding to m^'^ — 150 GeV), ti -[- Wt and 
14^14^ +jet events (overlayed, from back to front). 
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150 

mn [GeV] 



ee 


niH 


=150 GeV ma 


=180 GeV niH 


=200 GeV 


efi 


Signal 

Bkg: tt + Wt 
WW+iets 


+5.7 
-4 
-1 


-0.9 

+14 

+3 


-5.5 
+16 

+5 


MM 


Signal 

Bkg: tt + Wt 
WW+iets 


+2.6 
-2 



+0.5 

+11 




-3.4 
+16 
-1 


ee 


Signal 

Bkg: tt + Wt 
WW+iets 


+7.6 
-1 



-1.8 

+5 
+2 


-13.2 

+5 
+2 



TABLE II: Estimate of systematic shift 5m h [GeV], as defined 
in eq. ([6]). (Corresponding systematic error is different; see 
text.) 



m,H [GeV] 



e/i mode [GeV] 
^^ mode [GeV] 
ee mode [GeV] 



Combined [GeV] 



150 



4.1 




<3.1 



"STB" 



180 




< 8.2 



200 




< 11 



TABLE III: Estimates of the statistical error Amn'"*' in the 
niH reconstruction using Ni leptons. (e, fj, from t are not 
included.) The factors in the square-root correspond to unity 
for an integrated luminosity of 100 fb~^ and using leptons of 
the signal events. (iVf' = iV,'^"' + iV,'""' + N'f"\) 



FIG. 4: lO"" X (Cg) as a function of niH in the theoretical 
prediction for dV h^Uvij / dEo\Mii:<i5GcM ■ Muon energy distri- 
bution given in Fig. [3] is used: solid (dashed) line represents 
the expectation value taken with respect to the signal-plus- 
background (signal) MC events. 



close to 2Mw pJl], the correct solution can be identified 
relatively easily, using profiles of the lepton energy distri- 
bution, dilepton invariant mass distribution, etc. Hence, 
we consider deviations only around the correct solution. 
Fig. [3] shows the muon energy distributions separately 
for the signal events and the ti+ Wt and VFVF+jet back- 
ground events, which passed through all the cuts for the 
H — ^ jjLjjLUU mode. The signal events are generated with 
m^*-' = 150 GeV and the number of events after the 
cuts is about 12,000. The normalizations of the back- 
ground events have been rescaled according to the re- 
spective cross sections to match the number of the sig- 
nal events. Using this distribution the expectation value 
(Og) is computed; see Fig. |3l[l^ The observable Oq 



in eq. ([3]) is defined with G = G^^l^ and the theoreti- 
cal prediction for dTH-^ieuu/dEo is computed imposing a 
cut Mu < 45 GeV on the dilepton invariant mass. The 
value of (Og) changes as a function of mH-, which enters 
as an external parameter in the theoretical prediction for 
dT H^Ui^i^ / dEo\i\i^i,^4^QcY ■ Without any cuts and with- 
out background contributions, (Og) would cross zero at 
niH — 150 GeV in the large statistics limit. |20| As seen 
in Fig. m due to the effects of the cuts, the expectation 
value (Og) taken with respect to the signal events crosses 
zero at about -1-2.6 GeV above the MC input value. This 
value is consistent with Smn determined from eq. ([S]). 
{Og} with respect to the signal-plus-background events 



-0.8 GeV above the MC input 



crosses zero at Sttih 
value. dH 

In a similar manner, we can estimate the systematic 
shift SniH, as defined in eq. (jHl), in the rriH reconstruc- 
tion by varying the input Higgs mass value in MC simu- 
lations and using different modes. The values of Stuh for 
some sample points are listed in Tab. [Ill The magnitude 
of Smn due to all the cuts on the signal events is smaller 
for ££ = niJL than for e/i and is the largest for ee. This is 
because the acceptance corrections are smaller for /x than 
e. We confirm \5mH /'rnH\ ^ 1, which shows that devi- 
ations from the ideal limit is suppressed and kept under 
control. 

Tab. lIIII shows estimates of the statistical uncertainties 
(standard deviations) in the ran determination, Ato|^^'', 
in the case that Ni leptons are used to take the expec- 
tation value {Og)- The factors in the square-root corre- 
spond to unity for an integrated luminosity of 100 fb~^ 
and using leptons of the signal events in the respec- 
tive modes. The bottom row lists the upper bounds 
of combined statistical errors of the three modes; these 
are the largest among the three modes, when the num- 
ber of leptons is set as the sum of all three modes, 
^tot ^ Nl^'^^N^'^^+N'f^l We use Og with G = G^^a 
in these estimates. The estimates are derived in the fol- 
lowing way. We compute 50% = {0%) — {Og)'^ and 
d{OG)/dmH using the MC signal events which passed 
all the cuts. Since the statistical error of {Og) is given 
by AOg'^*- = SOg/V^I, we convert it to Amg'^*- using 
the tangent d{OG)/dmH evaluated at each input value 
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FIG. 5: 10^ X {Og) vs. mn (parameter in dV u^Uw I dEf\Mti<'&GcN^- MC events with tti^ = 150 GeV and for /i/.i mode are 
used. Through (a)-(e) only the signal events are used, while in (f) both signal and background events are used. All the cut 
values, except the ones shown explicitly in figs. (b)(c)(d)(f), are kept fixed to their default values. 
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Now we examine stability and reliability of our predic- 
tion, using the /i// mode and MC input m!^ — 150 GeV. 
In Fig. [5l^a) we vary the factorization scale /Xfac in PDF 
and PYTHIA and plot (Og) as a function of ran: /Xfac 
is taken as 1/2, 1 and 2 times its default value of MadE- 
vent (Pt of each scattered parton from each proton in the 
vector-boson fusion case).|22| Sran (and hence the value 
of nxH where (C'G('7iff))=0) changes by about 0.8 GeV 
as ^fac is varied from 1/2 to 2. In Fig. Eljb) we vary the 
cut value of the invariant mass of tagged two jets from 
550 GeV to 1 TeV. Corresponding variation of bTUH is 
about 0.8 GeV. In Fig. [SJc) we vary the values of the 



that (Og) as a function of vtih is also fairly stable. 

On the other hand, brriH and {OcijnH)) depend 
strongly on the lepton cuts. In Fig. [Sfd) we vary the val- 
ues of the Pt cut of the leading muon between 10 GeV 
and 30 GeV (while keeping the Pt cut value of the sub- 
leading muon to 10 GeV). The line {00(171^)) moves 
downwards and changes shape as the cut becomes tighter. 
This stems from the fact that the Pt cut suppresses lower 
part of the lepton energy distribution. 

Cuts involving missing transverse momentum P™"'''(= 
-Pt) may also affect the niH determination, since 



phad 



40 GeV + 6P and ^^^ *"^^^ affect the lepton energy distribution through in- 



Pt cuts of tagged jets, Pt ji > 

Pt,j2 > 20 GeV + 6PT,j, between 6Pt,j = and 2d'GeV. direct restrictions on PJjf. Here, PJjf denotes the trans- 



{SPt.j = is the default value.) Variation of Smn is 
about 1.0 GeV. These features agree with our expecta- 
tion that niH determination is insensitive to PDF and 
cuts involving only jets, since only the /3 distribution of 
the Higgs boson would be affected. Moreover, we find 



verse momentum of the tagged dilepton system, while 
P^^'^ denotes the sum of the transverse momenta of all 
other visible particles. Since magnitudes of systematic 
uncertainties in P^'^'^ and P^^ measurements are rather 
different, instead of varying cuts involving P™'*'*', we mul- 



tiply jPJi^'^l by a scale factor of 0.9, 1.0 and 1.1; we keep 
all the cut values unchanged. See Fig. [SJe). This varia- 
tion of energy scale affects restrictions on P^^ indirectly, 
since bounds on P^^^^ = —P^^'^ — P^^ are kept fixed. We 
find that variation of Smn is about 1.6 GeV. 

In Fig. \^i) we include the background contributions 
and vary the leading lepton Pt cut. (Through Fig.[5fa)- 
(e) only leptons from the signal events are used.) Com- 
pared to Fig. [SJd), each line moves upwards and the 
line shape is modified slightly by inclusion of the back- 
grounds. We have already seen this effect in the shift of 
the reconstructed Higgs mass in Fig. 2] and Tab. [TTl 

From these examinations and similar examinations for 
other input Higgs mass values and decay modes, we esti- 
mate that our prediction is fairly insensitive to uncertain- 
ties in PDF and jet variables. [23] On the other hand, the 
prediction is strongly dependent on the lepton Pt cuts. 
It is dependent also on other lepton acceptance correc- 
tions. These effects with respect to only leptons can in 
principle be estimated accurately, by understanding de- 
tector coverage and detector performances well. 

We note that all the lines in Fig. [Sjf) can be plot- 
ted using the real experimental data and can be com- 
pared with the prediction. Since these lines can be pre- 
dicted accurately and are dependent on MC input Higgs 
mass value nv^ , we can determine the Higgs mass by a 
fit of (00(771/^)), provided the background contributions 
can be estimated accurately. Similarly, all the lines in 
Figs. [SJb) and (c) can be compared, after inclusion of 
background contributions, with the corresponding ones 
plotted using the real experimental data. Alternatively 
we may make similar comparisons by loosening various 
cuts. This procedure raises relative weight of background 
contributions, so that we can test the prediction for the 
background contributions. We may also make use of the 
large degree of freedom of the observable Og for further 
tests of the prediction. 

In the case that higher lepton Pt cuts are unavoidable, 
determination of the Higgs mass by a fit of {OaimH)) 
would be more realistic, as mentioned above, since 
{Oci'mH)) is shifted considerably. [2J| Explicitly, a fit can 
be carried out in the following way. We can compute 
{OcimH)) using lepton energy distributions measured 
in a real experiment and MC simulation, respectively, as 

(OG(mH))exp = I dE,D,^^{Et,m'^^'')OG{Ei-m,H), (7) 

(Og("ih))mc = j dEiDMc{Ei;iv^'^)OG{Et,mH). (8) 

Defining a distance between {OG{'mH))cxp and 
(Og("ih))mc as 

/•ni2 

d^{m^'^)= / dmHw{mH) 

'J mi 

X [(Og(?«h)>cxp - {OG{mH))MC? ,(9) 

and minimizing d{m^'~^), we may reconstruct the Higgs 
mass. Here, w{mH) > denotes an appropriate weight. 
As the lepton Pt cuts become tighter, sensitivity to 



the Higgs mass decreases. For instance, in the case 
m%'"' = 150 GeV, if we impose Pr.iepi > 25 GeV and 
Pt, icp2 > 10 GeV to the leading and subleading lep- 
tons, respectively, we estimate that the statistical error 
of the Higgs mass, reconstructed by the above fit, to be 
about 5.0 GeV for an integrated luminosity of 100 fb~^. 
The size of error is not very dependent on choices of toi 
and 1712 or the weight w. (We took mi = 120 GeV, 
TO2 = 240 GeV and w = 1.) 

Comparing Tab. IIIII and the examination in Figs. [5l 
we anticipate that the statistical error will dominate 
over systematic errors. The accuracy of mn determina- 
tion is better for a lighter Higgs mass, within the range 
150 GeV <mH< 200 GeV. (e.g. If mn = 150 GeV, 
it would be reconstructed with 2-3% accuracy with 
100 fb~^ integrated luminosity.) 

Let us comment on the effects of leptons from r's, 
which we have neglected. Since these leptons would have 
a lower energy spectrum, the effects of the lepton Pt cuts 
are likely to be enhanced. We anticipate from Fig. [S{d) 
that the effects are to move the line (C'G("T-ff)) down- 
wards. Nevertheless, we expect that the additional con- 
tribution from t's and the cut effect will be accurately 
predictable and should not deteriorate the mn determi- 
nation. 

The main purpose of the present study on mn recon- 
struction is to demonstrate that the mass reconstruction 
using Og is applicable to such a complicated process. 
We chose the W^M^-fusion mode (rather than gg-insion 
mode) since a good signal-to-noise ratio is beneficial in 
demonstrating clearly the characteristics of the present 
method. We confirmed stability of our prediction ex- 
cept against purely leptonic cuts. Good understanding 
of background contributions is also important. 

Experimentally mn would be reconstructed very ac- 
curately using the decay modes H -^ ZZ'-*'> — > ££i£ in 
the relevant mass range [l3|. The lyVF'*-' decay modes 
can be used to test consistency of the Higg s mass val- 
ues, using our method or those of [ll|, U^- Ref. [ll[ 
uses leptonic variables but depends severely on PDF. We 
reemphasize that our method is independent of PDF in 
the leading-order. In a future work, we may examine the 
gg-fusion mode, in which the statistical error would be 
smaller while systematic errors would be larger than the 
T^T^-fusion mode. 

The observables proposed in this paper would have 
wide applications. For example, in supersymmetric mod- 
els, Og can be applied to a slepton decay into a lepton 
plus an undetected particle; since it is a two-body decay, 
its analysis would be straightforward. Another possible 
application is a measurement of the top quark mass. A 
preliminary analysis indicates that the top quark polar- 
ization effects are small. 

The authors express their condolences for victims of 
disasters in Japan. The disasters also struck the current 
project heavily. The authors hope the completion of this 
work to be a step forward in revival from the tragedy. 
The work of Sumino is supported in part by Grant-in- Aid 
for scientific research No. 23540281 from MEXT, Japan. 
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y is defined with respect to the boost direction of X. This 
should be distinguished from the (pseudo-) rapidity rj, de- 
fined with respect to the beam direction of experiments, 
used in the simulation study of mn reconstruction. [24] 

We note that if the shape of the lepton energy distribu- 
tion is unchanged, i.e. 5D oc D, 5D does not affect the 
reconstructed mass value. 

m_r/-dependent cut on the transverse mass Mt of U- 
P™'"" system (used in Tab. 7 of '6]) is omitted. We ex- 
pect many additional ways for optimization in the mn 



determination, such as inclusion of this cut. 
After this paper was submitted to arxiv, the ATLAS and 
CMS collaborations announced new Higgs mass bounds 
at the Lepton-Photon 2011 Conference [3|. Although the 
parameters ttih = 150, 180, 200 GeV are excluded at 
95 % C.L., the results are still preliminary and the exclu- 
sions are not conclusive. In particular, rriH ~ 150 GeV is 
at the boundary of the 95 % C.L. excluded region and it 
would be premature to exclude this case. 
The statistical error of reconstructed mn increases grad- 
ually with n. 

This mass region is being excluded by the recent Tevatron 
searches 8] ■ 

In Figs. |4]and[5] we choose the normalization of Oc as 
AT = l/r{H -^ fifiw; M^^ < 45 GeV) . 
The statistical error corresponding to the current MC 
events (A^^ « 24, 000) is Amg^*' ~ 0.6 GeV; c.f. Tab.llHl 
In this p articular example, however, since Srri'^ ^ 



-A%-f«if we can only estimate Smp+'°''^ ~ 0.8±0.6 GeV. 
The present statistics of the MC events should be suffi- 
cient for most of the other estimates of Smn in Tab. II. 
Our MC event generator being leading order, there is a 
significant scale dependence, for instance, in the distri- 
bution of the number of jets. 

This feature persists with higher lepton Pt cuts. For ex- 
ample, with a cut Pr.iepi > 25 GeV, Figs.[5la)(b)(c)(e) 
look qualitatively similar, except that all the lines move 
downwards. 

Recently requirement for the single lepton trigger has 
become severer in the LHC experiments 0]. Hence, pre- 
dictions with higher lepton Pt cut values may be more 
realistic. In principle, we can also use the dilepton trigger, 
in which case lower Pt cut values may be used. 



